On the simplified tree graphs in gravity 
A. I. Nikishov * 



P.N.Lebedev Physical Institute, Moscow, Russia 
November 4, 2011 

Abstract 

Firstly, I give the reason why is wrong my previously made assumption that the 
volume integral over the pressure may not be zero in a system where the gravitation 
plays no role in holding the system together. Secondly, in the first nonlinear approxi- 
mation I obtain the inner and outer Schwarzschild solutions in harmonic and isotropic 
coordinates in two different ways. One way is to start from standard solution and make 
the appropriate coordinate transformation. The other way is to use the perturbation 
theory with elements of Schwinger and Weinberg source approach. This latter method 
is applicable in general case and it is useful to study all its peculiarities on known sim- 
ple example such as Schwarzschild solution. It turns out that this method is simpler 
then S-metrics approach (previously made by Duff) and more informative as it it shows 

which contribution comes from what region of space. 
*E-mail: nikishov@lpi.ru 



1 introduction 

If the Earth gravitational field can be measured with G 2 — accuracy, than its changes, due 
to building up the pressure, leading to earthquakes, can be monitored. For this reason 
calculating the Earth gravitational field with this accuracy is of interest. Thus simplifying 
the perturbation methods is essential. As a first step in this direction we reproduce the 
corresponding approximation of the Schwarzschild solution by simple means. 

In three-graviton vertex all gravitons enter symmetrically. It can be simplified in partic- 
ular cases. For two on-shell gravitons this has been done in [1]. In our case one graviton is 
external field graviton. Using integration by parts we can release it from being differentiated. 
Then it interact with Weiberg gravitational energy-momentum tensor, see §6 Ch. 7 in [2]. 
In static case this tensor gives (negative) gravitational energy density which is three times 
that of Newtonian one [3]. Yet, as shown in [3], for a pressure-free dust there is an over all 
agreement with Newtonian limit for total energy of the system. We shall see later on that 
the same is true for the ball of liquid. 

The negative energy density makes me uneasy. I prefer to think that this is only effective 
quantity. In Newtonian theory the gravitational energy density is given either by or by 
— , see §99 in [4]. I think that the correct expression should have the form fJ,(f)+ in 
which the gravitation field has positive energy density. The relativistic generalization of this 
is given in [3]. 

In Section 2 I give the reason why is wrong my previously made assumption that the 
volume integral over the pressure may not be zero in a system where the gravitation plays 
no role in holding the system together. 

In Section 3 the interior Schwarzschild solutions in harmonic and isotropic coordinates (in 



G 2 — approximation) are obtained from the standard solution. They have the forms different 
from those obtained by using tree graphs or solving Einstein equations by perturbation 
procedure. Yet these differences are only in forms. 

In Section 4 the Weiberg gravitational energy-momentum tensor and relevant expressions 
for uniform liquid of a ball are given. 

In Section 5 the nonlinear corrections are obtained. 



2 Gravitational field of a system in which gravitation 
plays no role 

For a statical system the metric at large distance (in comparison with its size) is determined 
by the integral over energy-momentum tensor. In particular we need the integral f a^dV 
where is the stress tensor. This integral is equal to f PiXkdf where Pi is the force acting 
on the surfice of the body and the integral is over the surfice, see §3 in [5]. As in our case 
Pi = 0, the integral is zero and pressure plays no role in generating the gravitational field at 
large distances This assertion correct the mistaken assumption in [6]. 



3 Interior Schwarzschild solutions in harmonic and 
isotropic coordinates 

Here we obtain the solution in harmonic coordinates from the known solution in standard 
coordinates. The metric have the form, see Synge [7], eq.(7.183) 

r 2 dr 2 2 2 . 2 2 3y/l - qa 2 - y/l - qr 2 2 2 2mG 
ds =- + r z (d9 z + sin' 9dip) - {— — } dt , q= — — . (1) 



1 — qr 2 2 a 



In this Section we denote the harmonic coordinate radius by R and use the substitution 



mC R 2 

r = R(l-<l>(m.b.R)), <f>(m.b.R) = ^(-3 + (2) 

The radius of the liquid ball a in standard coordinates is related to that in harmonic coordi- 
nates b by the expression 

a = b + mG. (3) 
We are working in G 2 — approximation. It is easy to obtain 

dr 2 . mG (n R 2 , m 2 G 2 A 3 R 2 7R\ lln2 „,„ 3 . 
1 - qr 2 b b 2 b 2 9 2 6 2 4 b 4 

In terms oc G 2 we have replaced a by b. Using the relations 

RdR = X a dX ai (dR) 2 = XaXp ^ adXfi , a, (3 = 1,2,3., (5) 

we find, see (2) and (3) 

/ -rt\ , _,,n P Tfl 2 G 2 ,3R 2 2R 4 .X a Xa . . 

g afi (X) = [1 - 0(m, 6, i?)] 2 ^ + -^-(-^ yr ] ^^- (6) 

For cjr o we obtain from (1) 

r- „ . / \ ?n 2 G 2 .3 3 i? 2 3i? 4 N1 
Soo = -[1 + 20(m,a,r) + _(- - - — + -— )]. (7) 

Here again we have replaced a by b and r by i? in terms oc G 2 . Expressing a and r in terms 
of b and R we find 

, mG, o r\ , „, 3m 2 G 2 m 2 G 2 R 4 

20 m, a, r = -3 + - = 20 m, 6, i? + — . 8 

a a 2 b 2 o° 

So, from (7) and (8) we have 

, , N m 2 G 2 d5 3R 2 R\^ 
g 00 = -{1 + 20(m, b, R) + -^-( T " ^ " ^ ( 9 ) 



In isotropic coordinate system in the considered approximation b remains the same. In- 
deed, if we denote (only in the formula below) the radius of the ball in isotropic system as b, 
then 

a = 6(1 + ^)2 ~b + mG, (10) 

g 00 in isotropic coordinates has the same form (9) and g a p can be obtained from (6) by adding 
pure gauge functions with appropriate coefficients. The gauge function has the form 

In what follows we replace R by r. We need three such functions 

(x a ) p = 6 a/3 , (x a r 2 )^ = 5 afS r 2 + 2x a xp, (x a r 4 )^ = 5 aP r A + Ar 2 x a x p , (11) 
Adding to r.h.s. of (6) the function 

Kr = ^-{\*<# - ^(<W 2 + + 26*(^ 4 + 4r 2 ^)}, (12) 

we get in isotropic coordinate system 

, r . mG, r 2 N m 2 G 2 Ah 3r 2 3r\. , , , 

9«p(rn,r)=5 a p[l + — (3--) + — (-- — + — )}, r < b. (13) 

In both harmonic and isotropic systems the metric and their first derivatives are continuous 
across the boundary of the ball, see also [9] and [11]. 

4 Weinberg gravitational energy-momentum tensor 

In general relativity the nonlinear source of metric is the Weinberg gravitational energy- 
momentum tensor, see §7 in [2] 

T ik = + to,, i, k — 0, 1, 2, 3, 



tik = ^c\-\ h ^ R{1) + \vikh lm RV + R m _ ^ ikR (% (14) 

Here Tik is energy -momentum tensor of matter, and are linear and quadratic in h 
parts of R, 

12 

g%k = Vik + h ik , h ik =hik + hik +0(G ) (15) 
% = diag(-l, 1,1,1) RW^Rp, R( 2 )=R i (2)i = r ? ik R i ( k 2) , 

R ik = ^jAk - h J i, jk ~ h J k ,ji + h ik ,/). (16) 

The expression for R \ k ' see eq.(7.6.15) in [2] (or in another form in [3]). 
Using the expression for t a p in (33) below, it is easy to verify that 

= -^q^^ = ^ «, /3 = 1, 2, 3. (17) 

Hence 

(T a/3 + t aP )p = p' a + {Mf)' a = 0. (18) 

This is the equation for equilibrium of a liquid ball, see §3 in [8]. For ji = const differentiation 
over x a gives V 2 p = —AnGfi 2 or in spherical coordinates 

Twice integrating over r from to r we find 

p(r) - p(0) = -^27rG//V. (20) 



As p(b) = we have from here 



p(0) = InG^b 2 , (21) 



so that 



^ = ^^-0 = ^=^(3-^). (22) 



Here we neglect the difference between b and a and difference between m and to, see below. 
The Einstein equations in this formalism have the form [2] 

Rfk ~ \^kR {1) = -8nG(T ik + t ik ). (23). 

To calculate U k we need hf k . It is the same in harmonic and isotropic coordinates: 

h ik = -2<f>5ik,, h=hi 1 = -40. (24) 

Here is the Newtonian potential, see (26) below. 
For the ball of liquid with uniform density we have 

Ti j = (pl + p^u 1 + p5i j , T ° = -fi9(b - r). (25) 

Here b is the radius of the ball in the considered coordinate system. Now we have 

m'G r 2 

0( m ',r) = — (-3 + -), r<b, (26a) 
m'G 

<j,{m\r) = — . r > 6. (266) 

Here 

m' = ^TTfib 3 . (27) 
We remind here that in general relativity the dressed mass is given by the expression [8], 

[9] 

m = ^Tifxa 3 , (28) 

where a is the ball radius in standard Schwarzschild coordinates. The relation a = b + mG is 
exact in harmonic coordinates and approximate in isotropic one. The bare mass is, see §100 
in [7], or [9] 

m = J dVfxe x ^ 2 = J dVn{\ + ^r 2 G) == to + (29) 



From (27), (28) and (3) we have [9] 



m =m(l — )=m (l rr — )• (30) 

00. 



Now, in terms of (see (26)) tensor £ ife in (14) has the form 
1 



- -^{^ifc[3(V0) 2 + 400 jQ J - 20 5 j0 jfc - 400 jifc - S(f)(f) jaa 5io5 ko }, (31) 



8ttG 



' V 2 = 47rG/i. 



The indices of (true in the sense of general relativity) tensor are raised and lowed with g 
not with rj [2]. So, from (19) we have 

T ° = -fi, Tj = 6 aP p, go0 T = ^ + 2^ = T^+Tg\ T ( 0) = /i, T^ = 2/^, (32) 

and from (25) 

*oo = -^(V0) 2 -6/i0, 
W = g^{^[3(V0) 2 + 400, 77 ] - 20 ja i/3 - 400, Q/3 }. (33) 
We also need the quantities 

ioo = ^oo - ^Voot = ^~^( V< ^) 2 ~~ M 

t = t aa - t 00 = ^(W) 2 + 1O//0. (34) 



Too = Too - ^ooT = ^ + + ii<t> = T ( 0) + T 00 - (35) 



i 



Using (26) for r < 6 we get 



1 m 2 G 2 , A 20r\ 1 m 2 G 2 , t 21r\ 



1 m 2 G 2 , r 



2 



coo 



8nG 6 4 v 6 2 



(36) 



t ^_L^ {M _ 12 + £ } (37) 

'» = ? + ^=lSG^-^ (39) 



Similarly, for r > b we have 



3 m 2 G 2 - 1 m 2 G 2 5 m 2 G 2 

to °-"8^G^' '""^"H"' '"^"P"' (40) 

1 m 2 G 2 

Q/3 ~ 8ttG r 4 1 a/? ^ ( ^ 

Knowing t o (see (36) and (40))and (see(38)) we can calculate the total gravitational 
energy U : 

f „ r , ™riK / 69 3 12 18 x m 2 G 3m 2 G , . 

f/ = y dy (tM + Tfi') = (- - 5 - - - -)— = --— (42) 

The first number (i.e. ||) comes from t o when r < b, the second (— |) when r > b, the third 
term comes from and the last one from ml, see eq.(30). In the Newtonian theory we 
have 

.1 1 m 2 G 3m 2 G 
U =-(-2 + ^— = "5— ' < 43) 

The number | comes from ,r > b and from r < b. In terms oc G 2 we neglect the difference 
between m and m . 



Nonlinear corrections to metric 



2 



Up to gauge terms hik (rn',r) can be obtained from the formula (cf. §17 in [10]) 



2 



hik K, x) = / dV£> + (x - x')^(x'), (44) 



i 



^(x') = l6nG(T ik (x') + ^(rr')), (45) 



Here 



D+(x) = —5+(x 2 ) = - X — - 1 . 
+ v ; 4tt ; (2n) 2 x 2 + ie' 



di&D + {x) = -S(x). 



From (44), (46) we have 



djd 3 hik (m',r) = -9 ik (x), 



or in our case when hik (m',r) is independent of time 



V 2 hik (m',r) = -8 ik (x), 



We note also that 



/ 



drD + (x — x', r) 



47r\x — X'\ ' 



see eq.(3.92), Ch.2 in [10]. 
Using also the relation 

f 1 dt 

J -i \Jr 2 — r' 2 — Irr't 

it is easy to find that for any function /(r) we have 



2 if r' < r 

r ' 

if r' > r, 



/ d 4 x'D + (x-x')f(r') = - dr'r' 2 f(r') + dr'r'f(r'). 

J r'>c T J c J r 



Now for r < b, adding t 00 in (36) and T o in (39) we have 



1 _- i 1 m 2 G 2 5r 2 , 

-0 O o - too+ T oo - Y^-^( 9 " 



16nG 



Using also relations 



J d A x'D + {x-x')6{b-r') 



\b 2 -\r 2 , iir<b 



— if r > b 



J d 4 x'D + (x - x')r' 2 6(b - r') 



\b*-l Q r\ if r<b 



t *r>b, 



2 

we find the contribution to hoo from r' < b 

f j4 .„ , ^G 2 , 5r' 2 /r m' 2 ^ 2 ,^ 3 r 2 r\ 

J d*x'D + {x - — (9 - -p)9(b - r>) = -^-( T -a^+ft*)' r < 6 ' 

We note that r.h. sides of (53) and (54) are continuous at r = b 

— /i \ _ 2 

For r' > b we have T 00 = and too is given in (40). So the contribution to hoo 

16ttG / (fx'D + (x - x')—— = —r^, r < b. 
J r '>b 47rGr /4 b z 

Addind (55) and (56), we find 

2 , . , m 2 ^ 2 ^! 3r 2 r\ 

hoo (m,r) = - + ^' r < k 

Next, we express ft, ^(m',r) in terms of m, see eqs. (24), (26) and (30) 

i i m 2 G 2 r 2 

froo (m', r) = -2<f)(m', r) =hoo (m, r) + — — (-9 + — ), r < b. 



Thus, 



1 2 1 2 

/ioo (m ,r)+ /ioo (ra ,r) =hoo (m,r)+ hoo {m,r), 



where in agreement with (9) 



2 . . m 2 G 2 , 15 3r 2 r 4 , 
h»oM = —{-- + — + — ), r<6. 



2 

For r > 6 the contribution to /i o from r' > b is 



16.G / ^(x - *<)^ = ra>GH± - 1), r > 6. 



The contribution from r' < 6 is 



dVl^a; - a;')^-(9 - ^V)0(& - r') = m 2 G 2 4, r > b. 
o 4 tr rb 



Thus, the sum of these two contributions is 

koo (rri, r) = m 2 G 2 {\ -\), r > b. (63) 
rb r z 

On the other hand, from (24), (26) and (30) we have 

i . , . , f . 2m'G 2mG / _ i 3mG. 

/ioo (rri, r) = -20 m', r) = = (1 — , r > b (64) 

r r o 

Adding (63 and (64), we see that eq.(59) holds also for r > b. 
From (59) and (60), we have 

x m 3 G 2 , 15 3r 2 r 4 , 
0oo(»n, r) = -1 - 20(m, r) + — — (-— + — + — ), r < 6. (65) 

This agrees with Rosen [11], when terms up to G 2 are retained in his eq.(15). 
For r > b, we find 

2mG 2m 2 G 2 , , , 

g 00 (m,r) = -H — — , r > 6. (66) 

2 

To calculate h a /3 ( m ' \ r ) for r < 6 we note that it should have the form 

2 / / x r m 2 G 2 , r 2 r 4 . 

h a p(m,r)=6 a p ^ (c + c 2 ^ + c 4 — ).. (67) 

From here with the help of (57) we find 

22 2 m 2 G 2 21 3 r 2 1 r 4 

/i (m',r) =haa (m',r)- hoo (rri ,r) = — — {(3c - — ) + (3c 2 + -)^ + (3c 4 - 4)^4}, (68) 

and using the relation 

(r n ), aa = (n 2 + n)r n - 2 (69) 

for n = 2, and n = 4 we get from (68) 

P 2 m 2 G 2 r/ 3 X 6 , l N 20r 2 , N 

V 2 h (rri, r) = -^{(3c 2 + -) - + (3c 4 - ^-^}- (70) 



From (67) we obtain 



2 77l 2 G 2 T 2 

(h at3 (m',r)) iCe p = — ^(6c 2 + ^c A —\ 



Thus, 

r W(m',r) = V 2 h (m',r) - (h a p K,r)), a/3 = ^^{(12c 2 + 9) + (40c 4 - 5)-}. 

2 

On the other hand, we can obtain R W from the Einstein equation 

R = gnGTj = 8wG(-fi + 3p). 

In our case 

R =R R A (2) - K k R H\ R {1) = -2<p, aa = -8nG». 

So, 

R {1) = 8vrG3p- A (2) + fc* i2 2>, R « = 800, aa = ^(-36 + ^), 
r (2) = _ 8#>Qa _ 1O ( V 0) 2 = -^(36 - — ), . 
Collecting all these terms, we obtain 



6 4 v 6 2 
Comparison with (72) gives c 2 = —6, C4 = 3/4. 

2 

Now, using (67) and (69), with these c 2 ,C4 we can find the source of h a p (m',r): 

^ 2 2 / / x <■ m 2 G 2 . nn 15r 2 
V hap {m ,r) = d a p &4 (-36 + -^-). 

2 

For r < b this source contributes to h a p {m',r): 

f lA lry , Ar m 2 ^ 2 , 15r' 2 /7 ,, r m 2 G 2 ,57 6r 2 3r\ 

y dV£> + (* - ^)^_(36 - -gr^b - r>) = <W-^H T ~-gr + W )- 



For r > b its contribution is, see (53) and (54) 



d x D + (x - x )8 a p—^—{36 — ) 6{b -r) = 95 al3 — , r 



br 



For r > b we have (see also [3]) 



t a p = t aP --rj aP t= -^G m G (— — )■ r>fe - 



But this source generates the function, see [3] 



h a p= m G (— ^— ), 



which differs from the isotropic one by a multiple of the gauge function 



A Q ,/3 — m G ( — —). 



From the relation 



V 2 A a ,, = m 2 G 2 (-2^ + ^) 



we see that adding to t a p in (81) the gauge source 

2 

gives the source of h a p (m',r) in isotropic coordinates: 

^ = -3m 2 G 2 ^, r>b. 

2 

Thus , for r' > b and r < b the contribution to h a 8 from this source is 



I d 4 x'D + (x - x ')m 2 G 2 ^^- = -35 a sm 2 G 2 f°° dr'^- = ~6 
J r >>b r4 Jb r4 2 



3. m 2 G 2 

a/3 ^ 2 , 



Finally, adding this to (79) gives 



2 m 2 G 2 , 51 6r 2 3r\ 

W (m ,r) = ( T - + ^4), r < 6. 



1 

Using (24), (26a) and first eq. in (30) we express h a /3 (rn',r) in terms of m: 
W (m , r) =h a p (m, r) + ^—^—(-9 + -p-J, r < 6. 



From this and (88) we have 

12 1 2 

W (m , r)+ hap (m , r) =h a p (m, r)+ h a p (m, r), r < b. 



Here 



2 . . r m 2 G 2 Ah 3r 2 3r\ 

W {m, r) = 5 a p—^{— - — + —), r<b. 

For r > b the contribution from the source (86) is, see (51) 

/ d 4 x'D + (x-x')m 2 G 2 ^l = ^m 2 G 2 {- [ r dr' r -^+ f° dr'^-J 

Jr'>b r r Jb r Jr r 



-35 a ,m 2 G 2 {^ - ±}, r>b. 

Adding (80) and (92) we find 

2 6 3 

ha/3 (rri, r) = 5 a pm 2 G 2 {— + — }, r > b. 

i 

Using the first eq. in (30), we may rewrite h a p (m',r) in (24), (26b) as follows 

i . . N . 2m'G . ,2mG „ 2 ^ 2 1 , 

ha/3 {m,r) = d a/3 = b af} { 6m G —), r > b. 

r r rb 



From the sum of the last two equations we have 



Here 



12 1 2 

ha/3 (m , r)+ h a( 3 (m , r) =h a/ 3 (m, r)+ h a p (m , r). r > b. 



2 . . r 3m 2 G 2 

ha/3 (m, r) = b al3 ^ , r > 6. 



6 Concluding remarks 

One specific feature of the considered approach is that the metrics at any point is formed by 
the nonlinear sources in all space. In standard approach the inner and outer Schwarzschild 
solutions are obtained separately in each region. 

The second feature is the appearance of mass m! which depends on coordinate system. It 
can be expressed via invariant "dressed" mass m and then the solutions take a more simple 
form [9]. 

The gravitational energy density too of the Weinberg tensor is not the Newtonian one. 
This may indicate that in the future the gravitational theory will be modified to make the 
gravitational energy density positive. 
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